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Abstract 

A random walk starts from the origin of a <i-dimensional lattice. The 
occupation number n(x, t) equals unity if after t steps site x has been 
visited by the walk, and zero otherwise. We study translationally in- 
variant sums M(t) of observables defined locally on the field of occupa- 
tion numbers. Examples are the number S(t) of visited sites; the area 
E(t) of the (appropriately defined) surface of the set of visited sites; 
and, in dimension d = 3, the Euler index of this surface. In d > 3, the 
averages M(t) all increase linearly with t as t — > oo. We show that in 
d = 3, to leading order in an asymptotic expansion in t, the deviations 
from average AM (t) = M (t) — M{t) are, up to a normalization, all 
identical to a single "universal" random variable. This result resem- 
bles an earlier one in dimension d = 2; we show that this universality 
breaks down for d > 3. 
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1 Introduction 



We consider the simple random walk on a (i-dimensional hypercubic lattice: 
The walk starts at time t = at the origin x = and steps at t = 1,2,3,... 
with equal probability to one of the 2d neighboring lattice sites. The set of 
the sites that have been visited after t steps is called the support at time t. 
It is actually the support of the field of occupation numbers n(x, t) defined 
by 

, . I if x has not yet been visited at time t ,„ 
n (x,t) = , , . 1.1 
II otherwise 

The main variable that characterizes the support is its total number of 
sites S(t). The probability distribution of this random variable is known 
explicitly only in dimension d = 1, although it is well characterized [Q] also 
in higher dimensions. Our interest here will be in the connection between the 
statistical properties of S(t) and those of other variables characterizing the 
support. An example of the variables of interest is the (d — l)-dimensional 
surface area E(t) of the support, defined as the total number of nearest- 
neighbor bonds connecting a site of the support with a site not in it. We 
may imagine to obtain E{t) by inspecting all pairs of neighboring sites on 
the lattice and retaining a contribution unity whenever a pattern is found 
that consists of a visited and an unvisited site. Expressions for S{t) and 
E(t) in terms of the occupation numbers are easily constructed. For later 
convenience we shall express them instead in terms of the complementary 
occupation numbers 

m(x,t) = 1 -rt(x,i) (1.2) 

This leads to 

S(t) = £[l-ro(x,i)] 

X 

E(t) = ^2 [m(x, t) + m(x + e, t) - 2m(x, t)m(x + e, t)] (1.3) 

x e 

where e runs through the d basis vectors of the lattice. 

Both S(t) and E(t) count the total number of occurrences in space of 
specific local patterns of visited and unvisited sites. The quantities inside the 
summation signs in Eq. ( |Q| ) are the indicator functions of these patterns. 
Clearly observables can be constructed that count the total number of local 
patterns of a great many other types. A pattern may consist, for example, 
of a "dimer" of two visited nearest-neighbor sites oriented along a specified 
axis; or of an unvisited site having 2d visited neighbors; etc. etc. We shall 
generically denote such observables and their linear combinations by the 
symbol M(t). The class of all M(t) is defined precisely in Sec.|2[ 

There certainly exist many nontrivial correlations among the random 
variables M(t). Nevertheless, it should be possible, in principle, to express 
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them in terms of independently fluctuating degrees of freedom of the sup- 
port. Examination Q of dimension d = 2 has yielded a most surprising 
result. It appears that to leading order as t — > oo the deviations from 
average AM(t) = M(t) - M(t) satisfy 

AM (t) ~ b M tilogt)-^ 1 ' 1 r}(t) (d = 2) (1.4) 

where r)(t) is a random variable of zero average and unit variance, bM a 
proportionality constant, and k m a nonnegative integer called the order of 
M (one has for example k$ = and kg = 1). Eq. ( |1.4| ) says that in the large 
t limit the whole class of random variables M(t) reduces to just a single 
fluctuating degree of freedom n(t). 

In the present work we ask how the above conclusion is affected by the 
spatial dimensionality d. We are not aware of any simple heuristic argument 



that leads to (L4) and that would provide an indication whether or not a 
similar reduction of degrees of freedom still holds in higher dimension. It 
should be expected that as d goes up, independently fluctuating degrees of 
freedom will appear. 

We show by explicit calculation that in dimension d = 3, to leading order 
as t — > oo, the fluctuating degrees of freedom still collapse to a single one in 
a manner very similar to d = 2. Explicitly, 

AM(t) ~ am (tlogi)^ (d = 3) (1.5) 

where is a random variable of zero average and unit variance, and om 
a proportionality constant. Eq. ( [1.5] ) holds in particular for the choice M = 
S. The fluctuations AS(t) have been thoroughly studied and we know H] 
therefore that the probability distribution of £(t) is Gaussian. The variables 
£(t) in d = 3 and r/(t) in d = 2 are universal in the sense that they do not 
depend upon M. Differences between d = 3 and d = 2 have to do, among 
other things, with the appearance of an "order" &m in d = 2 but not in 
d = 3. Also, rj(t) is known to be non-Gaussian. 

The way to derive the three-dimensional result is closely analogous to 
what was done in two dimensions. We calculate the correlation of two ob- 
servables M(t) and M'{t) and find 



AM(t)AM>(t) ~a M a M >tlogt + 0(U) (f = 3) (1.6) 

where au and om' are known proportionality constants. We then exploit 
the fact that the prefactor in Eq. ([O]) factorizes into an M and an M' 
dependent part. It follows that the normalized random variables 

£ w (i) = a^itlogty^AMit) (1.7) 



have £ M (t) = and £m(*)£m'(*) = 1 for all M and M'. Therefore the 
averages and variances of all differences £m (t) — £m'(*) vanish and the £m(*) 
must all be identical to one and the same random variable that we shall call 
£(i); whence (O). 
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We show, again by explicit calculation, that in all dimensions d > 3 the 
correlation between two arbitrary observables M and M' has the form 

AM(t)AM'(t) ~ CMM't (d>3) (1.8) 

with a prefactor cmm' that no longer factorizes. Hence for d > 3 the fluc- 
tuations of the M(i) cannot be described any more by a single random 
variable. 

Our principal result, described above, is of a general nature. In Sec.|| 
we use our calculational framework to study three examples of fluctuating 
observables M(t). One of these is the Euler index x(t) of the surface of the 
support of the three-dimensional random walk. We recall that the Euler 
index of a closed two-dimensional surface embedded in three-dimensional 
space (the surface of a "Swiss cheese") is 

X = 2(1- H + C) (1.9) 

where H is the number of handles of the surface and C the number of cavities 
enclosed by it. In Sec.||, for appropriately defined "handles" and "cavities", 
we construct an explicit expression for x(t) and show that it is a member of 
the class of M(t). We then determine its asymptotic time dependence and 
its correlation with other observables. 



2 General observables M(t) 

The general property or "observable" M(t) that we shall deal with can now 
be defined as follows. Let A be a finite (typically small) subset of lattice 
vectors. Then we write 

M(t)=J2»AM A (t) (2.1) 

A 

where the ha are numerical coefficients, where the sum is through all subsets 
A not equivalent by a translation, and where 

M A (t) = Y,m x+A (t) (2.2) 

X 

with 

m x+A {t) = H m(x + a,t) (2.3) 

Here is a convenient notation for the set of lattice sites {x + a | a 6 A}. 

When A = the RHS of Eq. pjl) should be assigned the value unity and 



Eq. ( 4.11 ) shows that M$(i) is the total number of lattice sites. As an 
example one may check that S(t) has ha = ±1 for A = and A = {0}, 
respectively, and ha = otherwise; and that E(t) has [iA = 2d, —2 for 
A = {0}, {0,e}, respectively, and HA = otherwise. We note that the 
expansion coefficients \xa should satisfy the relation J2a A*A = if the sum 



on x in Eq. (|2.2|) is to have a finite value in the infinite lattice limit. 
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3 Averages of observables 



3.1 Averages M(t) 

The expression (^^) averaged on all random walk trajectories represents the 
probability that the set x + A has not yet been visited by the walk. It is not 
surprising, therefore, that the average M(t) of interest to us is connected to 
the distribution of first passage times for that set. More precisely, we can 
express it in terms of the probability / X +A( r ) that the walker's first visit to 
the set x + A takes place at time r. Elementary manipulations lead to the 
explicit relation 

-M(z) = -T^~ E ^aF a (z) (3.1) 
1 Z 



in which 



F A (z) = J2k+A(z) (3.2) 



and where, here and below, any time dependent quantity X(t) that occurs 
will be associated with a generating function X(z) = J2uLo z t X(t). The 
first passage time probabilities can be expressed in the random walk Green 
function by standard methods. After some algebra [|j one finds 

P A(z) = j^ 7) G A \z) (3.3) 

in which 

G~ A \z) = E (G^)aa' (3-4) 
a,a'eA 

where G A is the matrix of elements G(a — a',z) with a and a' restricted 
to the set A, and G(x, t) is the probability that a walker starting at the 
origin occupies site x after t steps. For the later purpose of discussing the 
small- (1 — z) expansion it will be convenient to split G(x, z) up according to 

G(x,z)=G(0,z)-«7(x,z) (3.5) 

With the aid of a little algebra we find the analogous splitting 

G A (z) = G(0,z)-g A (z) (3.6) 

in which 

a,a'eA 

where gyi is the matrix of elements g(a — a', z) with a and a' restricted to 
the set A. Using Eqs. ([O]) and ( |3.3| ) in Eq. ( p.l[ ) and inverting that equation 
we obtain the final result for the average of a general observable M(t), 
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where the integral runs counterclockwise around the origin. The coefficient 
Ga(1) appears in the literature || as the electrostatic capacity of the set A. 



3.2 Averages M(t)M'(t) 



The calculation of the correlation AM(t)AM'(t) requires that we now con- 
sider the average of the product M(t)M'(t). The latter can be rewritten 
as 

M{t)M\t) = E VAv' B M Au{r+B) (t) (3.9) 

r A,B 

The observable inside the sum on r is of the same type as M{t) and M'(t), 
even if it is parametrically dependent on r. The average of (^^) can therefore 
be evaluated with the aid of Eq. (3J3). After separating the terms with A = 
or B = from the others we find 



r<+1 ^ A,B#0 r (3.10) 



in which 



Iab(*, z) = -= 7-r - , r - . r (3.11) 

GAU(r+B)(Z) G A {z) G B {z) 

This expression involves GAu(r+B)( z ), a quantity whose r dependence needs 
to be rendered more explicit if the sum on r in Eq. ( |3.10| ) is to be performed. 
This we are able to achieve only in the scaling limit r — > oo and z — > 1 
with ^ 2 = 2dr 2 (l — z) fixed. In that limit a calculation analogous to that of 
Ref. [01 for the two-dimensional case leads to 



G A (z) + G B (z)-2G(r,z) 



(3.12) 



Gau(t+b)(z) G A (z)G B (z)-G 2 (r,z) 
where 

G(r, z) ~ d d /\- d /\l - ^) rf/2 - 1 e 1 " d/2 ^/2-i(0 (3.13) 



The expressions of this section up through Eq. ( 3.1 1| ) are exact for all finite 



times t. Only in Eqs. ( 3.12j ) and ( 3.13| ) has the scaling limit been taken, 



is an infinite lattice implied, and does the dimensionality of space appear 
explicitly. 



4 Long time behavior 

4.1 Averages M(t) 

The long time behavior of M(t) is determined by the behavior of M(z) in 
the vicinity of z = 1. We need the small- (1 — z) expansion of Ga(z), hence of 
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(5(0, z) and g& (z). The expansion of G(0, z) is well-known j|, ||] and reads 



f G -|^(l-z) 1 / 2 + O(l-z) (d = 3) 

G(0,z) = <^ Go-^(l-z)ln T l I + C4 (l-z) + ((l-z) 3 / 2 ) (d = 4) 
I G -c d (l-z) + O((l-zf/ 2 ) (d>5) 

(4.1) 

where we have abbreviated Go = G(0, 1) and where the with d = 4, 5, . . . 
are positive constants. Since g(x, z) is a solution of 

E [f( x + z ) - 2 9 ( x > 2 ) + 5(x - e M , 2f) = -5 X) o (1 - ^)G(0, z 



2d 



(4.2) 



one may conclude that g A {z) defined in Eq. ( |3.7D has an expansion of the 
form 



g A (z) = g A {\) + (1 - z)g { j ] + o(l - z) 



(4.3) 



as 2 — ► 1, where g A is an j4 dependent numerical constant. The (1 
expansion of M{z) therefore reads 



z - 



M(z) = < 



- z)- 2 + ff^m 2 (l - z)-*/ 2 + 0((1 - z)- 1 ) 
mi(l - z)~ 2 + ^m 2 (l - z)- 1 log ^ 

+ (m 2 - c 4 m 2 1) )(l - z)- 1 + o((l - z)" 1/2 ) 
k mi(l - z)~ 2 + (c d m 2 + mW)(l - z)- 1 + 0((1 - z)" 1 ^) 



(4.4) 

in dimensions d = 3, d = 4, and d > 5, respectively; here the coefficients 



(i) 

ni '"2 



m„ , m 



are determined by the observable M through the relations 



in 



niM] = ~ E T^ry 



m 



(i) 



[M] = - E 



(i) 
VA9 A 

i=0 



(4.5) 



When no confusion can arise we shall not indicate their M dependence 
explicitly. For transforming back to the time domain we need the basic 
expansions, for s > 0, 



dz 



t*- 1 r 



2vri 7 z t+1 (1 - z) s T(s) 



1 + 



(4.6) 



dz 1 
27ri J (1 - . 



In 



1 - z T(s) 



log*- ^i(s) + O 



logt^ 



(4.7) 
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where ip is the psi function. The first integral is easily found by writing 
the integrand as a power series in z and the second one follows from it by 
differentiation with respect to s. The result is 

_ ( ™iH^m 2 i 1/2 + 0(l) (d = 3) 

M(t) = I mit + ^ I m 2 \nt + m 1 + ^m 2 - c 4 m 2 + rr$ + o^ 1 / 2 ) (d = 4) 
[ mit + mi + c d m 2 + + 0(i~ 1/2 ) (d > 5) 

(4.8) 

where 7 denotes Euler's constant. In Sec. || some of the coefficients appearing 
in Eq. ([4.8D will be calculated for specific examples of observables M{t). 



4.2 Fluctuations and correlations 

In order to study the fluctuations and correlations of two observables M(t) 
and M'(t) in the large time limit, we now have to evaluate the average 
( |3.10| ) in the limit t —* 00. The principal problem is to find the expansion 
of J2r Iab(?, z) as z — ► 1. Whereas Ga(z) and Gb(z) remain finite in that 
limit, Eq. (|3.13| ) shows that in all d > 2 the quantity G(r, z) vanishes when 
z — > 1 (with £ fixed). This suggests that different leading powers of 1 — z 
are associated with different powers of G(r, z). We therefore express Iab in 
the form 

I AB (r, z) = h(z)G(r, z) + ■ ■ ■ + I n ^{z)G n ~ l {v , z) + /<g(r, z) 

(4.9) 

in which we take for Ik(z) the coefficient of (5 fc (r, z) in a power series ex- 
pansion of the RHS of Eq. ( |3.11 ) in which for l/GAu(r+B)( z ) the RHS of 



Eq. (|3.12j ) has been substituted. Explicitly, 

h(z) = -2/G A (z)G B {z) 

h{z) = (G A (z) + G B (z))/G 2 A (z)G 2 B (z) 

Hz) = -2/G A (z)G 2 B (z) (4.10) 

Even though Eq. ( p.!2| ) holds only in the scaling limit, Eq. ( |4.9| ) is an exact 
identity if it is taken as the definition of the remainder f AB - It is now 
possible to find the small 1 — z expansion of ^ T I AB {r,z) in the following 
way. After substituting Eq. ( [4.9[ ) in Eq. ( p.9| ) we evaluate the sum on all 
space of the powers G k (r,z) for k = 1,2,.... For k = 1 conservation of 
probability trivially leads to J2 r G(r,z) = (1 — z)~ l . For k = 2 one easily 
derives the general relation 

J2G 2 (r,z) = (l + z±)G(0,z) (4.11) 
after which use can be made of Eq. (|4.1| ). It follows that 
TG 2 (r z) = { |£(l--)- 1/2 + 0(l) 
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in dimensions d = 3 and d = 4, respectively. In dimensions d > 4 this 
quantity does not diverge as z — ► 1. Finally, for k = 3, we employ the scaling 
expression ( 3.13j ) for G(r,z). In dimension <i = 3 the resulting integral 



on £ diverges logarithmically at the origin, and using the lower end cutoff 
£ ~ ry/1 — z we find 

£G 3 (r,z) = ^ In ^ + 0(1) (d = 3) (4.13) 

In dimensions d > 3 the sum with fc = 3 remains finite as z — > 1. The re- 
mainder of the work will be carried out separately for the three qualitatively 
different cases d = 3, d = 4, and d > 5. 

4.2.1 Dimension d = 3 

We employ Eq. ( |4.9|) with n = 4. In the limit z — ► 1 the sum ^ r /^(r, 1) 
remains finite. Using the above expansions and also expanding GU^) and 
Gb(z) for small 1 — z we find for the average of Eq. fl3,10| ) the expression 



1 r dz 



M(t)M'(t) = — i — - 2m 1 m' 1 



1 



3 



2vri 7 L x 1 (1 - z ) 
3 5/2 / / 1 

+ nq/9 ( m l m 2 + Tl 1 m2) 



23/27r V-i'» 2 , -l"^ (1 _ z)5 /2 

+ £ m2m/2 7r 2 ^ log r 1 i + ---] (414) 



in which the are determined by M' via Eq. flj.5|) . After carrying out 
the z integral with the aid of the basic relations of Eqs. ( f4.6|) and (|4.7|) and 
subtracting the product M(t)M'(t) obtainable from Eq. (|4.8| ) we are led to 
the final result 

3 3 



AM{t)AM'(t) ~ _m 2 m£ tlnt + CMM't (4.15) 

where the constant cmm' is due to contributions from various subleading 
terms that have not been displayed explicitly above. Eq. ( |4.15 ) leads to the 



result announced in Eq. (|Oj) of the Introduction, whose consequence is the 
universality embodied by Eq. (|1.5|). 



4.2.2 Dimensions d = 4 and d > 5 

We employ Eq. (|4.9f ) with n = 3. In the limit z — > 1 the sum ^ r f^ B (r, 1) 
remains finite. Proceeding as in the case of d = 3 we obtain in dimension 
d = 4 the final result 



AM(t)AM'(t) = 

mim'i — Go(mi?7i2 + m2m' l ) + mim^' + mi — Fmm 1 



* + "(4.16) 
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where 

Tmm*= E MaMb £/S(r,l) ( 4 - 17 ) 

The only difference that arises when d > 5 is that not only J2 r /abC 1 ") 1) but 
also E r G 2 (r,l) remains finite as z — > 1. Straightforward calculation leads 
in a slightly different way to exactly the same result ( 4.16 ). If we denote the 
coefficient of the term linear in t again by cmm 1 -, then Eq. ( [4.16 ) constitutes 
the result announced in Eq. (|1,§|). 



5 Examples 

We apply the preceding theory to three examples. As a preliminary we ren- 
der more precise the concept of "surface of the support". We imagine M. 
partitioned into unit cubes that surround the lattice sites. Visited (unvis- 
ited) unit cubes correspond to visited (unvisited) lattice sites. The support 
of the walk then defines a visited volume whose (d — l)-dimensional surface 
will be called the surface of the support. In the second example we shall 
study the Euler index of the surface of the support of the three-dimensional 
walk. In order for this topological invariant - which involves the handles 
and the cavities of the support - to be well defined, we adopt the convention 
that two visited unit cubes are connected as soon as they share a vertex (and 
hence a fortiori when they share an edge or a face) . Fig. 1 shows an example 
of two connected unit cubes; their connectedness may be made visible by an 
infinitesimal deformation of the surface near the common vertex. 



5.1 Surface area E(t) of the support 

From Eq. (JfH^) we deduce that the average surface area of the support in 
dimension d > 3 behaves as 

E(t) = mi[E]t + o(t) (5.1) 

In the particular case d = 3 we have, more precisely, 

E(t) = mi[E]t + TJ ^m 2 [E}Vt + 0(l) (d = 3) (5.2) 



From Eq. ( |4.5| ) and the remarks below Eq. Q2.3| ) it follows that for all d 

m n [E]=2d[(G -~)- n -Gv n ] (5.3) 

where we used that for A = {0}, {0,e} one has Ga = Go, Go — \, respec- 
tively. In d = 3 we have G = 1.5164... (f§) so that 

mi[£] = 1.95, m 2 [E] = 3.20 (d = 3) (5.4) 

In dimension d 3> 1, using the expansion Go = 1 + (2<i) -1 + 0(d~ 2 ) one sees 
that E(t) ~ [2d - 3 + ©(cT 1 )]^ which may also be found by more heuristic 
arguments. 
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5.2 Euler index \{t) 

The Euler index of an arbitrary, not necessarily closed, two-dimensional 
surface embedded in M (with d > 2) is 

X = 2(l- J ff + C)-B (5.5) 

where -ff, C, and S are the number of handles, cavities (connected compo- 
nents of the complement of the surface), and boundaries, respectively. 

The support of the two-dimensional walk is a flat two-dimensional surface 
with H = C = and therefore has Euler index \ = 2 — B; the number of 
boundaries B is equal to one plus the number of unvisited islands enclosed 
by the support. This quantity was studied in Ref. ||. 

Here we wish to study the Euler index of the surface of the support of 
the t/iree-dimensional walk. For this closed (B = 0) two-dimensional surface 
the Euler index reduces to x = 2(1 — H — C). The handles and cavities are 
nonlocal objects of varying sizes and shapes, and it is interesting that one 
should be able to calculate the Euler index. The reason that this is possible 
is that the Euler index can be expressed, by means of the Gauss-Bonnet 
theorem, as an integral over the surface of a function of the local curvature. 

We first need a discrete analog of the Gauss-Bonnet theorem, that is, 
an expression for x as t ne sum over space of a function defined locally 
on a two-dimensional surface built out of plaquettes. In order to find the 
contribution from the surface around a specific vertex, we examine all 2 8 
patterns of occupation numbers in the surrounding cubic volume of 2 x 2 x 2 
unit cubes. The local contribution to the Euler index will be a weighted 
sum of the indicator functions of these patterns. There are only 22 patterns 
not equivalent under symmetry operations, and therefore 22 weights to be 
determined. One can set up equations for these weights by requiring that the 
final expression correctly yield the known result for various surfaces having 
the topology of a sphere (x = 2) or of a torus (x = 0). One finds, for 
example, that the pattern shown in Fig. 1 has weight —3. 

It is then a tedious but straightforward matter to cast the expression for 
X in the form of Eq. (|2.l| ) by working out the products with factors 1 — m, 
summing on x, and rearranging terms. Since the resulting expression for 
X may also be of use in other contexts, we have listed the coefficients in 
Fig. 2. The nonzero coefficients ha are the 18 entries of that table, and 
those obtainable from them by rotations and reflections. Also listed are 
the values of the capacities Gvi(l), found numerically by determining the 
required values of G(x, 1) and inverting the 18 matrices G^. These data 
then allow the large time behavior of the average %(i) to be calculated. It 
has the general form of Eq. (4.8), 



2 l/2 3 3/2 

X(t)=m 1 [ x ]t+ m 2 [ X ]Vi + 0(l) (d = 3) (5.6) 



with the numerical values 

mx[x] = -3.44 x 10" 2 , m 2 [x\ = 2-12 x KT 1 (d = 3) 



(5.7) 
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The negativity of mi[x] implies that the Euler index is dominated by the 
number of handles. 

The Euler index is defined also for (d — l)-dimensional surfaces and can 
always be expressed as the integral of a local function on this surface. We 
have found no general way to obtain discrete analogs of these expressions 
when the surface dimension is d — 1, with d > 4, even though for a given 
dimension this is in principle possible. 



5.3 Cavities 

A quantity of potential interest is the number of cavities of a given size, shape 
and orientation. The indicator function of such a cavity is a single product 
of m's and (1 — m)'s and their total number is obtained by summing that 
product on space. This number therefore belongs to the class of obervables 
M(t) and behaves for large times as shown in Eq. (|4.8| ). 

Let C\(t) be the number of cavities consisting of a single site. We deter- 
mined the numerical values 

mi[Ci] = 7.92 x 10 -4 , m 2 [C 1 ] = -7.10 x 1(T 3 (d = 3) 

(5.8) 

These show the rarity of such cavities: Their number increases at the ap- 
proximate rate of one every one thousand steps. Cavities larger than a single 
site will certainly be still rarer, so that C\{t) ~ C(t). It follows that x(t) 
is approximately equal to minus twice the number H{t) of handles of the 
surface of the support. 



5.4 Correlations 

We have not included in the above examples the number of sites S(t) of 
the support, as this quantity has been thoroughly studied in the literature 
|], It is now useful to observe that the average S(t) behaves, in d = 3, 
according to the general law (4.8) with coefficients mi [5] and m^S] that 



are easily determined via our general rules. One finds 

m x [S] = 0.659, m 2 [S]= 0.434 (d = 3) (5.9) 

in agreement with the literature [f|j . Upon inserting the various coefficients 
77i2 determined above in Eq. ( 4.1 5| ) we directly obtain the leading order be- 



havior of all cross-correlations between the quantities studied. It appears 
that the total surface area E(t) of the support is positively correlated with 
S(t), and that (since C\(t) « C(t)) the numbers of handles and of cavities, 
H(t) and C(t), are negatively correlated with S(t). 



6 Discussion 



The asymptotic behavior as t — > oo of the average and the variance of the 
total number S of sites in the support is well known in all dimensions. 
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These known results are reproduced by the relations of this paper if one sets 
M = M' = S. In the dimensions d = 2 and d = 3 the variable S is known 

1/2 

to have anomalous fluctuations in the sense that the rms deviation AS 2 
increases faster with t than proportional to S (t). The fluctuations in S 
do not arise, therefore, as the sum of effectively independent contributions 
along the random walk trajectory, but must be due to a collective effect. It is 
easy to accept that similar anomalous fluctuations appear in an entire class 
of observables M(t) related to the support, and in their correlations. What 
is surprising, however, is that the M(t) all fluctuate in strict proportionality 
to one another. This proportionality holds both in dimension d = 3 (a result 
of this work) and in d = 2 (derived in Ref. 0). In d > 4 the anomalous 
fluctuations and the rigorous proportionality of the fluctuating observables 
disappear. 

The number of sites S(t) visited by a d-dimensional lattice random walk 
has its analog, in continuum space, in the d-dimensional volume (" Wiener 
sausage") swept out during a time interval t by a hypersphere of radius 
b that performs Brownian motion. The volume of the Wiener sausage has 
been a recurrent object of investigation (see e.g. Ref. ||]), and the islands en- 
closed by it in d = 2 have been considered by mathematicians {e.g. recently 
by Werner H). However, we are not aware of any results, in dimensions 
d > 3, on the total surface area of the Wiener sausage, let alone on such 
properties as the number of handles or cavities. We see here problems for 
future research. 
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FIGURE CAPTIONS 



Figure 1 . Example of a local configuration of 2 x 2 x 2 sites that may occur 
in the support of the three-dimensional walk. Shaded (unshaded) unit cubes 
are centered on visited (unvisited) lattice sites. The volumes of shaded cubes 
are considered to be connected via the common vertex. 



Figure 2. The Euler index can be written as x = J2a HaMa, with the Ma 
defined in the text. The first and fifth column list the subsets A of lattice 
sites (heavy dots) that have nonzero ijla- The other nonzero /ia correspond 
to the sets A obtainable from those shown by rotations and/or reflections. 
The total number of distinct sets related to A by such symmetry operations 
is OA- The coefficients Ga(X) are the capacities needed in the calculation. 
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